Met ast ability in zero-temperature dynamics: 
Statistics of attractors 

C Godrechet and J M LuckJ 

f Service de Physique de I'Etat Condense, CEA Saclay, 91191 Gif-sur-Yvette cedex, 
France 

I Service de Physique Theorique|, CEA Saclay, 91191 Gif-sur-Yvette cedex, France 

Abstract. The zero-temperature dynamics of simple models such as Ising ferro- 
magnets provides, as an alternative to the mean-field situation, interesting examples 
of dynamical systems with many attractors (absorbing configurations, blocked confi- 
gurations, zero-temperature metastable states). After a brief review of metastability in 
the mean-field ferromagnet and of the droplet picture, we focus our attention onto zero- 
temperature single-spin-fiip dynamics of ferromagnetic Ising models. The situations 
leading to metastability are characterized. The statistics and the spatial structure of 
the attractors thus obtained are investigated, and put in perspective with uniform 
a priori ensembles. We review the vast amount of exact results available in one 
dimension, and present original results on the square and honeycomb lattices. 



PACS numbers: 05.70.Ln, 64.10. +h, 64.60.My 

E-mail: godrecheOdrecajn. saclay. cea. f r,luck@spht . saclay. cea.fr 



t URA 2306 of CNRS 



Metastability in zero-temperature dynamics 



2 



1. Introduction 

The non-equilibrium dynamics observed in a variety of systems, ranging from glasses 
to granular media, with its characteristic features of glassiness and aging [T1I2], is often 
thought of as the motion of a particle in a complex energy landscape, with many valleys 
and barriers p|. This picture fully applies in the mean-field situation, where barriers 
are infinitely high and valleys appear as metastable states [H |3] . For finite-dimensional 
systems with short-range interactions, barrier heights and lifetimes are always finite at 
finite temperature, so that metastability becomes a matter of time scales jHlIZ]. 

The main goal of the present paper is to emphasize that the zero-temperature 
dynamics of simple models such as Ising ferromagnets provides, as an alternative to the 
mean-field situation, interesting examples of dynamical systems with many attractors 
(absorbing configurations, blocked configurations, zero-temperature metastable states). 
For completeness, we first give a brief review of the phenomenon of metastability 
in the mean-field ferromagnet (Section |21), and of the droplet picture implying that 
barrier heights become finite for short-range interactions (Section El). We then turn 
to an overview of the configurational entropy in disordered and complex systems, 
and of the so-called Edwards hypothesis (Section Ej). The rest of this paper is 
devoted to zero-temperature single-spin-fiip dynamics of ferromagnetic Ising models. We 
characterize the dynamics which lead to zero-temperature metastability, and investigate 
the structure of the attractors thus obtained. After a general presentation (Sectional), 
we review the extensive amount of exact results available in one dimension ([Hlin], and 
especially ^UI) (Section IH)), and present original results on the square and honeycomb 
lattices (Section [7j). We end up with a brief discussion (Section IS]). 

2. Metastability in the mean-field ferromagnet 

The concepts of metastability and of spinodal line date back to the early days of 
Thermodynamics, with the seminal work of Gibbs ^1], in the context of phase separation 
in simple systems such as fluids. In order to summarize these early developments, let 
us use the more modern language of the Landau theory for a ferromagnet At 
the mean- field level, i.e., neglecting any spatial dependence of the order parameter, a 
ferromagnetic system of Ising type is described by the Landau free energy density 



where T is the temperature, Tc the Curie temperature, A and B are positive phenome- 
nological parameters, and H is the applied magnetic field. The magnetization M{T, H) 
is such that F{M) is stationary. It is therefore given by the equation of state 




M2 



(2.1) 



—— = A(T - T^M + BM^ -H = 0. 
The corresponding magnetic susceptibility x is such that 




(2.2) 




(2.3) 
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At fixed temperature T <Tc, the solutions to (|2.2|) are of tlie following three types: 

(I) stable: absolute (global) minimum of the free energy, for \M\ > Mq{T), 

(II) metastable: relative (local) minimum of the free energy, for Msp(T) < |M| < Mo(T), 

(III) unstable: maximum of the free energy, for |M| < Msp(T). 





M 




H 


• 

• 




-M (T) 


• 

• 


T-T 

c 

M 


/ m 

/ 1 

1 1 

-M„(T) -M^^(T) 


\\ I 

m \°\ 

\ \ 
\ \ 



Figure 1. The three types of solutions to the equation of state (|2.2j) . Upper 
panel: Plot of the magnetization M{H, T) against the magnetic field for 
a fixed temperature T < T^- Thin full line: stable solution (I). Thick full 
line: metastable solution (II). Thick dotted line: unstable solution (III). Lower 
panel: Phase diagram in the M~T plane. Phases I to III correspond to the 
three types of solutions. Full line: spontaneous magnetizations M = ±Mo{T). 
Dashed line: spinodal magnetizations M = ibMsp(r). 

Figure [T] illustrates this discussion. When the magnetic field H is positive, the 
stable solution M{H, T) describes the magnetization of the + phase. If H decreases to 
zero, the magnetization takes its spontaneous value M(0,T) = Mq{T). Coexistence is 
reached: the + phase becomes degenerate with the — phase. If H decreases further to a 
small negative value, the + phase still exists but becomes metastable. Its magnetization 
M{H, T) can be analytically continued all the way until the endpoint of the metastability 
region, which is reached for OF /DM = d'^F/dM'^ = 0. The corresponding values of 
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M and H are respectively called the spinodal magnetization Msp{T) and the coercive 
magnetic field Hco{T). All physical quantities are singular at the spinodal endpoint. 
The spontaneous magnetization, spinodal magnetization, and coercive field read 

Mo(T) = M(0,T)= (^^^^^^^ ^ , 

M.p(T) = M(-/7eo(T),r) = ( = (2.4) 



HUT) 



SB J y/3 



V 275 ; 3v^ 
The metastable solution is separated from the stable one by an extensive free 
energy barrier, whose height AF grows proportionally to the volume Q of the system. 
The activation energy density AF is the difference in F{M, T) between the metastable 
solution (II) and the unstable one (III). It has a finite limit 

at coexistence, whereas it vanishes at the spinodal endpoint. 

The continuation of the stable phase to the metastable one corresponds to an 
analytical continuation in the mathematical sense. The magnetization M{H, T), defined 
as the stable solution to ()2.2j) for > 0, is indeed analytic at H = 0. The power-series 
expansion 

M{H,T) = Y.an{T)H^ (2.6) 

n>0 

has a finite radius of convergence, equal to \Hco{T)\. 



3. Finite-dimensional ferromagnet: Metastability vs. droplet picture 

It has been realized in the 1960s that metastability is an artefact of the mean-field 
approximation. For realistic Ising ferromagnetic systems with short-range interactions 
in a finite-dimensional space. Fisher and Langer have shown that the 
magnetization M{H,T) of the + phase, or, equivalently, its free energy F{H,T), is not 
analytic at the coexistence point H = 0. Their argument relies on the droplet theory 
initiated by Frenkel ^Hl- Consider, in a (i-dimensional ferromagnetic model below its 
Curie temperature, a spherical droplet with radius R of the — phase (M = — Mo(T)) 
inside an infinite + phase (M = +Mo{T)). The corresponding excess free energy 

f{R) = 2ujR'^HMo{T) + dujR'^'^ao{T) (3.1) 

is the sum of a volume term, proportional to the volume uR"^ of the droplet, to the 
spontaneous magnetization Mq{T), and to the magnetic field H, and of a surface term, 
proportional to the surface area duR"^'^ of the droplet and to the interfacial free energy 
(To(r). Here u = 2%'^^'^ /T(d/2) denotes the volume of the d-dimensional unit sphere. If 
the magnetic field H is positive, the excess free energy is a positive and growing function 
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of the radius R. If the magnetic field is negative, however, the excess free energy reaches 
a maximum jT3] , 

/e = /(i?e) = ujaoirr ( ^J~\rr., ] , (3.2) 



2\H\Mo{T) 



for a finite critical radius 



_ (rf-l)ao(T) 

" 2|i7|Mo(T) • ^^-^^ 
The existence of a critical droplet with a finite radius Rc, and hence of a finite 
barrier height /c, has far-reaching consequences, both in statics and in dynamics. The 
reduced barrier height fc/T of the critical droplet (in units where Boltzmann's constant 
is unity) plays a role similar to that of the reduced instanton action S/h in Quantum 
Mechanics |16j . 

• Statics. The free energy F{H) of the + phase is not analytic at H = 0. In other 
words, the formal power-series expansion ()2.6|) is divergent: its radius of convergence 
vanishes. It can indeed be shown that F{H) has a branch cut starting aX H = 0, 
with an exponentially small imaginary part for if < 0, of the form 

Im F{H + iO) exp(^-j;^ . (3.4) 

• Dynamics. If the magnetic field is instantaneously turned from a positive to a small 
negative value, the 'metastable' + phase has a finite lifetime r, whose order of 
magnitude is the inverse nucleation rate of a critical droplet, given by an Arrhenius 
law in terms of the above barrier height: 

^ exp (-^] . (3.5) 



/ c 



Near coexistence {\H\ — ^ 0), the estimates ()3.4|) and ()3.5|) have an exponentially 
small essential singularity of the form 

i^,..(H.iO,^e.p(-^(^)"). ,3.) 

To sum up, for ferromagnets and similar systems, the metastability scenario stricto 
sensu (analytic free energy at the coexistence point, infinite lifetime in the metastable 
region, sharp spinodal line) is a peculiar feature of the mean-field and of the zero- 
temperature situations. For systems with short-range interactions, at finite temperature 
in finite dimension, the barrier height fc is finite all over the would-be metastable region. 
The free energy has an essential singularity at the coexistence point, and the lifetime of 
the metastable phase is finite, although it becomes exponentially large near coexistence. 
The spinodal line therefore becomes a crossover between the fast decay of an unstable 
phase (no barrier height) and the slow decay of an approximately metastable phase 
(finite barrier height). This crossover phenomenon has been recently investigated by 
means of a sophisticated numerical approach [T7] . 
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4. Metastable states and complexity in disordered and complex systems 

From the 1970s on, a variety of complex systems, such as structural glasses, spin glasses, 
and granular materials, have been shown to possess many metastable states at low 
enough temperature, or high enough density. The number Af{N; E) of metastable states 
with given energy density E typically grows exponentially with the system size, as 

Ar(iV;E)~exp(iV5eonf(^)). (4.1) 

The quantity Sconi{,E) is referred to as the configurational entropy, or complexity [TSj . 

Most investigations of the complexity were motivated by dynamical considerations. 
The dynamics of the above systems in their low-temperature or high-density regime often 
turns out to be so slow that the system falls out of equilibrium P, and exhibits aging 
and other characteristic features of glassy dynamics [2]. It has been proposed long ago 
to describe slow relaxational dynamics in terms of the motion of a particle in a complex 
energy (or free energy) landscape [3], with many valleys separated by barriers. Several 
approaches have been developed, in order to make this heuristic picture more precise, 
and chiefly to identify these valleys. Metastable states have thus been rediscovered 
under various names and with various definitions, in different contexts: valleys [S], pure 
states IHlin], inherent structures [201, quasi-states j2I]. 

In the mean-field situation, all these metastable states are separated by extensive 
barriers, so that they have an infinite lifetime, just as the unfavored phase in an 
Ising ferromagnet. For finite-dimensional systems with short-range interactions, barrier 
heights and lifetimes are always finite at finite temperature, so that metastability 
becomes a matter of time scales. Furthermore, from a dynamical viewpoint, the various 
concepts of metastability recalled above are not equivalent |S[7|. 

Whenever metastable states live forever, i.e., either in the mean-field situation or 
at zero temperature, a natural question concerns their dynamical weights. Consider 
a system instantaneously quenched into its low-temperature or high-density regime, 
from a randomly chosen initial configuration. Does the system sample all the possible 
metastable states with given energy density with equal statistical weights, i.e., with a 
uniform or flat measure, or, to the contrary, does any detailed feature, such as the shape 
or size of the attraction basin of each metastable state, matter? This question is of 
interest for any dynamical system having a large number of attractors. 

A similar question has been addressed in many works on the tapping of granular 
systems jHll22l2n]- Under tapping, a granular material constantly jumps from a blocked 
configuration to a nearby one. In this context, Edwards [24 proposed to describe 
situations such as slow compaction dynamics, or the steady-state dynamics under gentle 
tapping, by means of a flat ensemble average over the a priori ensemble oi all the blocked 
configurations of the grains with prescribed density. 

Extending the range of application of this idea far beyond its original scope, 
the so-called Edwards hypothesis commonly refers to the assumption that all the 
metastable states with given energy density are equivalent. This hypothesis has two 
consequences. First, the value of an observable can be obtained by a flat average over 
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the a priori ensemble, or Edwards ensemble, of all those metastable states. Second, 
the configurational temperature, or Edwards temperature, Tconf = (dS'conf/d-E)"^, is 
expected to share the usual thermodynamical properties of a temperature, and especially 
to coincide with the effective temperature involved in the generalized fluctuation- 
dissipation formula in the appropriate temporal regime [22] • The concept of ergodicity, 
and the resulting thermodynamical construction, therefore hold, as in equilibrium 
situations, up to the replacement of configurations by metastable states, and of 
temperature by the configurational temperature Tconf- 

The Edwards hypothesis has been shown to be valid for the slow relaxational 
dynamics of some mean-field models j2I], where metastable states are indeed explored 
with a fiat measure. This hypothesis has also been found to hold for finite-dimensional 
systems, at least as a good numerical approximation, both in the steady-state of tapped 
systems in the regime of weak tapping intensities [HI 122 EH], and in the regime of slow 
relaxational dynamics of various models [21]. 

Besides the mean-field geometry, another physical situation where metastable states 
are unambiguously defined is the zero-temperature limit, where no barrier can be crossed 
at all. The metastable states are then defined as the blocked configurations under 
the chosen dynamics. For instance, for an Ising model with single-spin-fiip dynamics, 
a metastable state is a configuration where each spin is aligned with its local field, 
whenever the latter is non-zero. The rest of this paper is devoted to the statistical 
analysis of the blocked configurations thus generated, before we come back to the so- 
called Edwards hypothesis in Section |H1 



5. Zero-temperature single-spin-fiip dynamics of the Ising model: 
Generalities 

This section is devoted to an overview of zero-temperature dynamics of ferromagnetic 
Ising models. In the two following sections, we review the extensive available results in 
one dimension ([Hlin], and especially JO]); and present novel results in two dimensions. 
The ferromagnetic Ising model is defined by the Hamiltonian 

^ = - XI ^mO-n, (5.1) 

{m,n) 

where the (T„ = ±1 are classical Ising spins sitting at the vertices of a regular lattice 
with coordination number z, and the sum runs over all pairs (m, n) of nearest neighbors. 

Consider single-spin-fiip dynamics (i.e., Glauber dynamics, in a broad sense) in 
continuous time. Each spin is flipped {(Xn — > — o"„) with a rate W{ST-C) per unit time. 
This rate is assumed to only depend on the energy difference implied in the flip, 

6n = 2hn(Tn, (5.2) 

where 

dT~C 

K = = '^rn, (5.3) 

^ m{n) 
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is the local field acting on spin n. The sum in the above definition runs over the z 
nearest neighbors m of site n. 

The condition of detailed balance at temperature T = 1/13 with respect to the 
Hamiltonian Ti reads 

= exp(-;3^W). (5.4) 

The two most usual choices of flipping rates, or flipping probabilities in the case of 
discrete updates, 

Heat-bath: Whb(5^) ^ 



l + exp(2/357^)' (5.5) 
Metropolis: Wui^H) = min(l, exp(-/357^)), 

obey the detailed balance condition ()5.4p at any temperature. It is, however, less well- 
known that there is a vast family of dynamical rates, besides these two choices, which 
obeys the condition ()5.4|1 . Indeed, the energy difference 571 may assume z+1 equidistant 
integer values: 

5n = 2z, 2^-4, -2z, (5.6) 

so that there are a priori z+1 different rates. If the coordination number z is even 
(resp. odd), the condition ()5.4|) only gives z/2 (resp. {z + l)/2) equations for the z + 1 
unknown rates, so that there remain [z + 2)/2 (resp. {z + l)/2) free parameters (one of 
them corresponding to the absolute time scale). 

The arbitrariness described above persists in the zero-temperature limit. The 
condition of detailed balance ()5.4|) indeed reads 

W{6n) = for all 6n > 0, (5.7) 

whereas nothing is imposed on the rates with STi. < 0. For the sake of deflniteness, we 
focus our attention onto the class of zero-temperature dynamics deflned by the following 
flipping rates: 

r if > 0, 

W{5n) = lWo ii6n = 0, (5.8) 
i 1 if6n<0. 

The rate Wo relative to the free spins is kept as a parameter. A free spin is deflned as 
a spin cr„ which is subjected to a zero local fleld (/i„ = 0, hence 6H = 0). Free spins 
can only exist if the coordination number z is even, so that belongs to the list ()5.6|) . 
The zero-temperature limits of the heat-bath and Metropolis rates ()5.5|1 are respectively 

VToHB = 1/2 and Wou = 1- 

The zero-temperature dynamics of ferromagnetic Ising systems is far from being a 
trivial problem in general 12H] ■ The single-spin-flip dynamics deflned by ()5.8|) is a 
descent dynamics (i.e., every move strictly lowers the total energy) in the following two 
situations: 



• z is odd, so that there are no free spins at all. 
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• z is even and Wq = 0. The corresponding dynamics is said to be constrained. As 
such it belongs to the class of kinetically constrained systems l3Uj . 

In both situations, every spin flips a finite number of times, and the system gets 
trapped in a non-trivial attractor (absorbing configuration, blocked configuration, zero- 
temperature metastable state). 



6. Constrained zero-temperature single-spin-fiip dynamics on the Ising 
chain 

The one-dimensional situation of a chain of spins corresponds to the smallest possible 
value of the coordination number [z = 2), so that the list of values of the energy 
difference is 4, 0, and —4. The heat-bath dynamics has been solved exactly by 
Glauber 

Let us focus our attention onto zero-temperature dynamics. For any non-zero value 
of the rate Wq corresponding to free spins, the dynamics belongs to the universality class 
of the zero-temperature Glauber model. This is a prototypical example of phase ordering 
by domain growth (coarsening) [H2] • The typical size of ordered domains of consecutive 
+ and — spins grows as L{t) ~ t^/^. The energy density E{t) = — 1 + 2/L{t) therefore 
relaxes to its ground-state value as The particular value Wq = corresponds to the 

constrained zero-temperature Glauber dynamics, which has been investigated in [31^], 
and especially in ^U]. Hereafter we summarize the main results of these investigations. 

6.1. Mapping onto the dimer RSA model 

In the constrained zero-temperature Glauber dynamics, the only possible moves are flips 
of isolated spins: 

-+ > , +-+ ^ + + +. (6.1) 

Each move suppresses two consecutive unsatisfied bonds. The system therefore 
eventually reaches a blocked configuration, where there is no isolated spin, i.e., up and 
down spins form clusters whose length is at least two. Equivalently, each unsatisfied 
bond (or domain wall) is isolated. These blocked configurations are the zero-temperature 
analogues of metastable states. 

The dynamics can be recast in the following illuminating way. Going to the dual 
lattice, where dual sites represent bonds, let us represent unsatisfied bonds as empty 
dual sites, and satisfied bonds as occupied dual sites: 

= a„a„+i = -l^o, ^^^^^ 

SO that the moves ()6.H) read 

o o — > • • . (6.3) 

This mapping shows at once that the dynamics is fully irreversible, in the sense that 
each spin flips at most once during the whole history of the sample. 
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The dynamics (j6.3p is identical to that of the random sequential adsorption (RSA) 
of lattice dimers, which has been considered long ago jHSllSll- This is one of the simplest 
examples of RSA problems, for which powerful analytical techniques are available in one 
dimension, both on the lattice and in the continuum [35J. A blocked configuration thus 
appears as a jammed state of the dimer RSA model, as illustrated on the following 
example: 



(6.4) 



6.2. Relaxation of the m,ean energy 



We consider as initial state the equilibrium state at finite temperature Tq. Each bond 
variable is independently drawn from the binary distribution 

J — 1 (o) with probability p, . 
" 1+1 (•) with probability 1 — 
where the parameter p is related to the energy density Eq of the initial state and to the 
corresponding temperature Tq by 

f5o = ^ = lln^, Eo = -l + 2p. (6.6) 
To 2 p 



It is a common feature of one-dimensional RSA problems ^5] that the densities 
of certain patterns obey closed rate equations. In the present case, the densities pe{t) 
of clusters made of exactly i consecutive unsatisfied bonds (empty dual sites) obey the 
linear equations [101 ESI 

^ = -{i - l)pe{t) + 2 Y: Pkit) (6.7) 

for f > 1, with p^(0) = (1 — pYp^. These equations can be solved by making the Ansatz 
Piit) = a(t) z(tY for i > 1. The solution thus obtained, 

p,{t) = (1 - pe-'f exp(2p(e-* - 1)) p^e'^^"^)*, (6.8) 

leads to 

E{t) = -l + 2j2^Peit) = -1 + 2pexp(2p(e-* - 1)). (6.9) 

£>1 

Only clusters of length £ = 1 survive in the blocked configurations, as could be expected. 
These defects occur with a finite density, 

pi(oo) =pe-2p, (6.10) 

so that the energy density of blocked configurations takes a non-trivial value, which 
depends continuously on the initial temperature via the parameter p ^ 

E^ = E{oo) = -l + 2pe'^P. (6.11) 

This non-trivial dependence is a clear evidence that the dynamics is not ergodic. For an 
initial state close to the ferromagnetic ground-state {Eq — > —1, i.e., p ^ 0), the behavior 
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Eoo ^ Eq — is easily explained in terms of the clusters made of two empty sites. 
The energy of blocked states then increases monotonically with p, up to the maximum 
value Eoo = — 1 + e~^ ^ —0.632121, for an uncorrelated (infinite-temperature) initial 
state (p = 1/2, i.e., Eq = 0), and then decreases monotonically with p, down to the 
value Eoo = —1 -|-2e~^ ~ —0.729329, corresponding to an antiferromagnetically ordered 
initial state {p = 1). 

6.3. Distribution of the blocking time 

The late stages of the dynamics are dominated by an exponentially small density of 
surviving clusters made of two empty sites. More precisely, ()6.8|) shows that their 
density reads P2{t) ~ ae"*, with a = j9^e~^^. The dynamics can therefore be effectively 
described by a collection of aN such clusters, each cluster decaying exponentially with 
unit rate, when a down spin flips. The blocking time T/v is the largest of the decay 
times of those clusters. For a large sample, it is therefore distributed according to 
extreme- value statistics [36J: 



where the fluctuation Xn remains of order unity, and is asymptotically distributed 
according to the Gumbel law: f{X) = exp(— X — e~"^). 

6.4. Distribution of the final energy: Dynamics vs. a priori ensemble 

The attractors of the zero-temperature constrained single-spin-fiip dynamics of the Ising 
chain are the spin configurations where every unsatisfied bond is isolated. The most 
natural statistical description of these attractors is provided by the a priori ensemble 
where all the blocked configurations with given energy density E are taken with equal 
weights. 

In this section, the exact distribution of the final energy of the blocked configura- 
tions is compared to the prediction of the a priori approach. It turns out that there are 
two a priori predictions, the second (refined) one being a far better approximation than 
the first (naive) one. 

For a finite chain made of spins, the number of blocked configurations with 
exactly M unsatisfied bonds, i.e., empty dual sites, reads 



where 6 = (resp. 1) for a closed (resp. open) chain. This is indeed the number of ways 
of inserting the M empty sites in the N — M + b spaces made available by the presence 
of — M occupied sites, with at most one empty site per available space [12]. In the 
limit of a large system (M, 00, with a fixed ratio M/N = [E + l)/2), this number 
grows exponentially, irrespective of b, as 



Tn = HaN) + Xn, 



(6.12) 




(6.13) 



Af{N; E) ^ exp{N S^piE)) 



(6.14) 



Metastahility in zero -temperature dynamics 12 
in accord with (|4.1|) . The a priori configurational entropy reads [HIEZIEHI 

S^^{E) = E ln(-2E) + hi(l -E)- ln(l + E). (6.15) 

This result can be alternatively derived by means of the transfer-matrix method ISHl 
HUj . For a finite chain of spins, we introduce the partition function 

c 

where the sum runs over all the blocked configurations C, as well as the partial partition 
functions Z^, labeled by the prescribed value of the last bond, according to ()6.2|) . 
The latter quantities obey the recursion 

-(e^ o)- 

The 2x2 transfer matrix T has eigenvalues 



e^±v/4+^ 

A± = ^ , (6.18) 

so that Z'j^ ~ Z'^ ~ exp(A^ln A+). The configurational entropy Ss_p{E) is therefore given 
by a Legendre transform: 

lnA+(7)+7i5^-Mi5^) = 0, ^ = -^^, ^=^' (^-^9) 

hence 

e^ —2E 

E = ^ 7 = In ^ =, (6.20) 

so that ()6.15j) is recovered. 

The most naive prediction for the probability distribution f{E) of the final energy 
density is that f{E) is proportional to the number A/'(A^; E). We thus obtain a large- 
deviation estimate of the form 

/(E)~exp(-iVSap(E)), (6.21) 

where the naive prediction for the large-deviation function Eap(-E') reads 

S,p(E)=5^ax-5ap(^), (6.22) 

with 

5max = In ^^^^ ^ 0.481212 (6.23) 

being the maximum of the configurational entropy Sap{E). This maximum is reached 
for 

^max = -4e ~ -0.447214, (6.24) 

which is therefore the typical a priori energy density of a blocked configuration. 

A more refined prediction for the distribution of the final energy density, introduced 
in the context of Kawasaki dynamics j3T] , resides on the following observation. The mean 



Metastability in zero-temperature dynamics 



13 



energy density of blocked configurations, as they are generated by the zero-temperature 
dynamics, is equal to E^q, given by ()6.11|) . which is in general different from the a 
priori value ()fj.24j) . Within the a priori formalism under consideration, this difference 
is taken into account by attributing to every configuration C an extra weight of the 
form exp(— 7effA^-E(C)). The effective inverse temperature 7efr is chosen so that the 
mean energy density coincides with ()6.11|) . hence 7eflf = 'j{Eao), where 'j{E) = dSa.p/dE 
is given in ()6.20|) . This procedure amounts to replacing the configurational entropy 
S'ap(-E), which is maximal aX E = -Emax, by the relative entropy 

S,p{E\E^) = S,^{E) -{E- E^)-f{E^), (6.25) 

which is maximal ai E = E^o by construction. The resulting prediction for the 
distribution of the final energy density now involves the refined large-deviation function 

S,p(E|Eoo) = ^ap(^oo) - 5ap(^) + (E - E^^E^). (6.26) 

It turns out that the distribution of the final energy density of the blocked 
configurations, as they are generated by the zero-temperature dynamics, has been 
evaluated exactly by analytical means JU] • We prefer to skip every detail of this rather 
lengthy derivation and to state the result. The distribution of the final energy density 
is given by an exponential estimate of the form 

/(E)~exp(-iVSdyn(i?)), (6.27) 

where the large-deviation function Sdyn(-E) depends on the parameter p characterizing 
the initial state, and reads, in parametric form: 

;i + {2p - l)zf -{z- 1) , l + {2p- l)z + (1 - z)e'^P' 



-'dyn 



^^^^ Apz'^e'^P' {2{1 - p) + {2p - l)z) 1 + (2p- 1)2 - (1 - ;z)e 



2pz '■ 



2pz^e^P'{2{l -p) + {2p - l)z) ^ ' 

Figure |21 shows a comparison between the dynamical entropy Sdyn(-E), given 
by ()6.28|) . for an uncorrelated (infinite-temperature) initial state {p = 1/2), and the 
predictions ()6.22j) of the naive a priori approach and ()6.26|1 of refined a priori approach. 
The refined prediction turns out to be a good approximation to the dynamical entropy, 
although it is not exact. The following particular values allow a quantitative comparison. 
At the ground-state energy density E = —1, one has 

Sdyn(-l) = Inze(p), S,p(-l|Eoo) = In ^ _ (6.29) 

with 

1 + (2p - l)zc + (1 - z^)e^P^^ = 0, (6.30) 

i.e., for p = 1/2, Sdy„(-1) ^ 0.245660, Sap(-l|^oo) ~ 0.255408. At the maximum 
allowed energy E = 0, one has 

1 1 1 - ve~'^P 

Sdyn(O) = -- lnp(l - p), S,p(0|Eoo) = - In (6.31) 
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Figure 2. Full line: plot of the dynamical entropy of the ferromagnetic chain 
with constrained Glauber dynamics, given by Ti^y^{E) l|6.28jl . against energy E, 
for an uncorrelated initial state {p = 1/2). Dash-dotted line: prediction H6.22() 
of the naive a priori approach. Dashed line: prediction (|6.26() of the refined a 
priori approach (after _10j. 



i.e., for p = 1/2, Sdy„(0) ^ 0.693147, T.^p{<d\E^) ^ 0.744940. 

Another quantity of interest is the dynamical specific heat Cdyn, which is defined 
as the curvature of the dynamical entropy at the mean energy density, 

Sdyn(i?) - (6.32) 

SO that the energy variance (in the sense of sample to sample fluctuations of the observed 
energy per spin in the blocked configurations) behaves as 

VarE^^. (6.33) 

The exact dynamical entropy (j6.28|) and the refined a priori prediction (j6.26|) respectively 
yield 

Cdyn = 4pe-^f (1 - p + 4p2), Cap = 4pe-2f (1 - p^'''^){\ - 2pe-^n, (6-34) 
i.e., for p = 1/2, Cdy„ ^ 0.406006, C^p ^ 0.379540. 



6.5. Correlations: Dynamics vs. a priori ensemble 

Let us define the spin correlation function C„ and the connected energy correlation 
function G„ as 

Cn = (o^oC^n), Gn = (toT^) - E'^ = (aoCriO-„0-„+i) - (6.35) 

This section is devoted to a comparison between the exact energy correlation 
function G„ and the prediction of the refined a priori description. It turns out that 
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the spin correlation function C„ is more difficult to handle, both from the a priori and 
from the dynamical viewpoint. 

The energy correlation function (Gn)ap in the a priori ensemble at fixed energy 
density E can again be evaluated by the transfer-matrix method. We have, for n > 
in the bulk of an infinitely long chain, 

(G„)ap = {L+\S\R-){L^\S\R+) (^^^ . (6.36) 

In this expression, S = diag(— 1, +1) is the energy operator, while 



(^±I = 3^(A± e^), |i?±)= (^"'^^^ ) (6.37) 



are the left and right eigenvectors of T associated with the eigenvalues X±. We have 
consistently {L+\S\R+) = E. After some algebra we obtain the following expression, 
involving only the mean energy E jHj: 

/ 1 -I- F \ " 

{GnU = {l-E')[-^) . (6.38) 

The connected energy correlation function in the a priori ensemble thus exhibits an 
exponential fall-off, modulated by an oscillating sign. 

The energy (i.e., occupancy, in the RSA language) correlation function of the 
blocked configurations, as they are generated by the zero-temperature dynamics, has 
also been evaluated exactly (TU]. Skipping again every detail, one has 

{G^W = 2pe-^^ f (1 - 2p)tMl. ^2p • (6-39) 

\ m>n+l J 

Whenever p ^ 1/2, i.e., Tq ^ oo, the ffist term is leading, hence (G'„)dyn ~ {—2pY/n\. 
In the case of an infinite initial temperature {p = 1/2), one has (G'n)dyn ~ (— l)"'/(n + l)!. 
The connected energy correlations therefore exhibit a factorial decay, modulated by an 
oscillating sign. 



7. Zero-temperature single-spin-flip dynamics in two dimensions: 
Honeycomb vs. square lattice 

Leaving aside the one- dimensional realm, where exact analytical results are available, 
by means of an exact mapping of the dynamics onto an RSA problem, let us now turn 
to some novel results pertaining to two-dimensional examples of the zero-temperature 
single-spin- flip dynamics described in Sectional We will investigate the following cases: 

• Honeycomb lattice. The coordination number z = 3 is odd, so that the standard 
zero-temperature dynamics ()5.8|) is a descent dynamics. 

• Square lattice. The coordination number z = 4 is even, so that the situation 
is qualitatively similar to that of the chain. The dynamics (j5.8p is a bona fide 
coarsening dynamics for any non-zero value of the rate Wq, whereas the constrained 
dynamics (Wq = 0) is a descent dynamics. 
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In the following we discuss in parallel the standard zero-temperature dynamics on 
the honeycomb lattice and the constrained one on the square lattice. Both dynamics 
lead to metastability: the system gets trapped in a finite time into one of many blocked 
configurations (zero-temperature metastable states). Figure El shows typical blocked 
configuration generated by both dynamics described above on samples of 150^ spins. 
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Figure 3. Typical blocked configuration generated by the zero-temperature 
single-spin-flip dynamics considered here. Black (resp. white) sites represent 
up (resp. down) spins. Left: dynamics on the honeycomb lattice (z = 3 is odd). 
Right: constrained dynamics on the square lattice (z = 4 is even, Wq = 0). 



7.1. Distribution of the number of flips 

The constrained single-spin-fiip dynamics on the Ising chain described in Section IHl is 
fully irreversible, in the sense that every spin fiips at most once during the whole history 
of the system, before a blocked configuration is reached. 

In the present situation, although the two dynamics considered are descent 
dynamics in a global sense (the total energy decreases at each spin fiip), a given spin 
can fiip more than once. It turns out that the maximum number of fiips is Ms = 2 
for the constrained dynamics on the square lattice, and Mh = 5 for the dynamics 
on the honeycomb lattice. Skipping details, let us mention that these bounds can be 
derived by considering the number of unsatisfied bonds around the four squares, or the 
three hexagons, which surround a given spin. Figure |3] shows an example of a spin 
whose history contains this maximum number of fiips in each case. The mean number 
of fiips of a given spin remains however well below unity, at least for uncorrelated 
initial configurations. For dynamics on the honeycomb lattice, the mean number of 
fiips is (M)h ~ 0.4244, while the fraction of spins which fiip five times is of order 
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Ph(5) ~ 10 ^. For constrained dynamics on the square lattice, the mean number of flips 
is (M)s ^ 0.2577, while the fraction of spins which flip twice is ps(2) ~ 0.00203. 

u^ti iMX^ iXt^ nxu riTn 

tmj mxt tuxj tmj trtt 

F F 




F F F 



F F 

Figure 4. Examples of histories such that a given spin do experiences the 
maximum allowed number of flips. At each step the spin about to flip is circled. 
The letter F marks those steps where ao is about to flip. Upper panel (left to 
right): square lattice: ao (at center) flips Ms = 2 times. Lower panel (top 
to bottom and left to right): honeycomb lattice: ctq (at bottom of uppermost 
hexagon) flips Mh = 5 times. 




7.2. Distribution of the blocking time 

We have measured by means of numerical simulations the blocking time of finite samples 
of size N X N for the descent dynamics on the honeycomb and square lattice described 
above, with uncorrelated random initial configurations. Figure El shows a plot of the 
mean blocking time (T/v) against the size (number of spins A^^) of the simulated samples. 
The data are well fitted by the form (T/v) = aln A^^ + b + c/N. 

This observation suggests that the late stages of the descent dynamics can again 
be described in an effective way by a dilute population of two-level systems, which relax 
independently of each other, with a common characteristic time r, so that their density 
falls off exponentially, as p{t) ~ exp(— t/r). Extreme- value statistics indeed implies 

{TN)^T\nN^. (7.1) 

The fits shown in Figure |31 yield the characteristic times th ~ 2.49 for the honeycomb 
lattice, and rs ~ 1.84 for the square lattice. The fact that these estimates do not 
identify with simple numbers, together with the rather large observed correction to the 
result (|7.1|) . fitted by a 1/N term, suggests however that there might be more than one 
type of two-level systems. 



Metastability in zero-temperature dynamics 18 

25 

20 

A 15 

Vio 

5 



"0 2 4 6 8 10 12 

In 

Figure 5. Plot of the mean blocking time (T/v) of zero-temperature single- 
spin-flip dynamics, against the number of spins N'^. Full symbols: standard 
dynamics on the honeycomb lattice. Empty symbols: constrained dynamics on 
the square lattice. Full lines: fits (see text) with asymptotic slopes th ~ 2.49 
and Ts PS 1.84. 

7.3. Spin correlations 

We have measured the on-axis spin correlation function C„ = (aoCTn), where n denotes 
the site at distance n lattice bonds from the origin along one of the main axes of either 
lattice, and where (...) again denotes an average over the blocked configurations. The 
spin correlation at unit distance Ci is nothing but minus the mean energy density E^o 
of the blocked configurations per bond. We obtain (Ci)h = — (-Eoo)h ~ 0.7485 for the 
honeycomb lattice, and (Ci)s = — (-E'oo)s ~ 0.6126 for the square lattice. 

Figure [Ul shows a logarithmic plot of the on-axis spin correlation function, against 
distance n. The data at large enough distances are well fitted by straight lines, 
demonstrating that correlations fall off exponentially to zero, as C„ ~ exp(— n/^). The 
inverse slopes of the least-square fits yield ~ 1.27 for the honeycomb lattice, and 
~ 0.79 for the square lattice. 

8. Discussion 

The main focus of the present paper is on zero-temperature single-spin-flip dynamics 
of ferromagnetic Ising models. We have emphasized that there are two kinds of such 
dynamics. The first kind (X-type dynamics, in the classification of [2H]) generically 
corresponds to bona fide coarsening dynamics, and lead to phase ordering by domain 
growth, whereas the second kind (JF-type dynamics) corresponds to descent dynamics, 
and provides interesting examples of dynamical systems with many attractors. 

The one-dimensional case has been reviewed in Section IHl The existence of an exact 
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Figure 6. Logarithmic plot of the on-axis spin correlation function Cn against 
distance n. Full symbols: standard dynamics on the honeycomb lattice. Empty 
symbols: constrained dynamics on the square lattice. Full straight lines: least- 
square fits with inverse slopes ~ 1-27 and ~ 0.79. 



mapping of the constrained zero-temperature dynamics onto the dimer RSA problem 
opens up the possibility of an exact analytical evaluation of many dynamical observables, 
including characteristics of the attractors jTUI. The comparison between the exact 
results thus obtained and the predictions of the a priori ensemble allows to test the 
validity of the so-called Edwards hypothesis, in a regime which is very distant from 
the situation of gentle tapping, considered originally by Edwards. On the one hand, 
a flat average over the a priori ensemble of attractors (especially in its refined version, 
defined by imposing the value of the observed mean energy density) provides inexact, 
but numerically reasonable, predictions for physical quantities. The numerical accuracy 
of a priori predictions may even by very good, like e.g. for the dynamical entropy (see 
Figure 121). On the other hand, the energy correlation function provides an example of a 
qualitative discrepancy between the actual behavior of attractors and the prediction of 
any a priori ensemble. Connected correlations indeed exhibit a factorial fall-off in l/n\ 
with distance, which is a generic characteristic of lattice RSA problems, whereas any flat 
measure over an a priori ensemble can be described by the transfer-matrix formalism, 
and therefore leads to exponentially decaying correlations in one dimension. Another 
example of a qualitative discrepancy concerns the lengths of ordered clusters in one- 
dimensional spin models. Cluster lengths are automatically statistically independent 
and exponentially distributed in a priori ensembles, again as a consequence of the 
underlying transfer-matrix formalism, whereas both of these properties have been shown 
to be violated in several examples of zero-temperature dynamics |38t I41j. 

Examples of zero-temperature dynamics on two-dimensional Ising models have 
then been investigated numerically. The full dynamics on the honeycomb lattice and 
the constrained one on the square lattice exhibit very similar features. The observed 
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logarithmic growth of the mean blocking time suggests that the late stages of both 
examples of descent dynamics are dominated by the relaxation of independent two- 
level systems, in analogy with the one-dimensional situation. This mechanism, in 
turn, is rather suggestive that attractors are likely to have a non-trivial statistics. 
The quantitative comparison with uniform measures on a priori ensembles is far more 
difficult than in one dimension, because a priori ensembles themselves may possess 
phase transitions as a function of their density. In particular, the accurately observed 
exponential fall-off of spin correlations is not very informative in that respect. Finally, it 
is worth noticing that the metastable states put forward in a recent work on competitive 
cluster growth j32] share many of the features of the present two-dimensional attractors. 

To sum up this investigation of the statistics of attractors of zero-temperature 
descent dynamics in finite dimensions, it turns out that the so-called Edwards hypothesis 
does not hold in this context in general. The predictions of the a priori approach, 
however, often provide good numerical approximations. This picture is corroborated by 
a recent work by Camia [IH] . 
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